Checkerboard density of states in strong-coupling superconductors 
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The Bogoliubov-de Gennes (BdG) equations are solved in tlie strong-coupling limit, where real- 
space (preformed) pairs bose-condense with finite center-of-mass momenta. There are two energy 
scales in this regime, a temperature independent incoherent gap Ap and a temperature dependent 
coherent gap Ac{T), modulated in real space. The single-particle density of states (DOS) reveals 
checkerboard modulations similar to the tunnelling DOS in cuprates. 
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Many independent observations show that the super- 
conducting state of cuprates is as anomalous as the nor- 
mal state. In particular, there is strong evidence for 
a d-like order parameter, which changes sign when the 
Cu02 plane is rotated by 7r/2 0. A few phase-sensitive 
experiments |3| provide unambiguous evidence in this di- 
rection. A d-wave BCS gap could appear in the two- 
dimensional Hubbard model near half filling, as suggested 
by Scalapino, Loh, and Hirsch j^] concurrently with the 
discovery of novel superconductors. On the other hand, 
c-axis Josephson tunnelling high- precision magnetic 
measurements ^ , photo-excited quasi-particle relaxation 
dynamics |^ and some other experiments support 
more conventional anisotropic s-like gap. 

In fact, there are stronger deviations from the conven- 
tional Fermi/BCS-liquid behaviour than the gap sym- 
metry. There is now convincing evidence for pairing 
of carriers well above Tc as predicted by the bipolaron 
theory Q, the clearest one is provided by the uniform 
magnetic susceptibility |^ 0|, tunneUing, and photoe- 
mission. The tunnelling and photoemissi on g ap is al- 
most temperature independent below Tc [Q ^^'^ 
exists above Tc pH Il3l IT^ with its maximum several 
times lar ger t han expected in the weak and intermediate- 
coupli ng JlSlI BCS theory. Kinetic 1171 and thermody- 
namic [13 data suggest that the gap opens both in charge 
and spin channels at any relevant temperature in a wide 
range of doping. A plausible explanation is that the nor- 
mal state (pseudo)gap, Ap, is half of the bipolaron bind- 
ing energv |18l. although alternative models have been 
proposed 

Further studies of the gap function revealed even more 
complicated physics. Reflection experiments, in which 
an incoming electron from the normal side of a nor- 
mal/superconducting contact is reflected as a hole along 
the same trajectory l2(|| , showed a much smaller gap edge 
than the bias at the tunnelling conductance maxima [21| . 
Two distinctly different gaps with different magnetic field 
and temperature dependence were observed in the c- 
axis I(V) characteristics They were also observed 

with the femtosecond time-resolved optical spectroscopy 
[2^ . More recent STM experiments revealed checker- 



board spatial modulations of the tunnelling DOS, with 
[2^ 1 and without j2^ I2HI applied magnetic fields. 

We have proposed a simple model |[27| explaining two 
different gaps in cuprates. The main assumption, sup- 
ported by a parameter-free estimate of the Fermi energy 
|28| . is that the attractive potential is large compared 
with the renormalised Fermi energy, so that the ground 
state is the Bose-Einstein condensate of tightly bound 
real-space pairs. In this letter I calculate the single parti- 
cle DOS of strong-coupling (bosonic) superconductors by 
solving the inhomogeneous BdG equations. When pairs 
are Bose-condcnscd with finite center-of-mass momenta, I 
obtain a checkerboard DOS reminiscent of the tunnelling 
DOS in cuprates. 

The anomalous Bogoliubov-Gor'kov average 



F,,,(ri,r2) = (#,(ri)*,,(r2 



is the superconducting order parameter both in the weak 
and strong-coupling regimes. It depends on the relative 
coordinate p = — r2 of two electrons (holes), described 
by field operators ^s(r), and on the center-of-mass coor- 
dinate R = (ri -|-r2)/2. Its Fourier transform, /(k, K), 
depends on the relative momentum k and on the center- 
of-mass momentum K. In the BCS theory K =0, and 
the Fourier transform of the order parameter is propor- 
tional to the gap in the quasi-particle excitation spec- 
trum, /(k, K) oc Ak. Hence the symmetry of the order 
parameter and the symmetry of the gap are the same in 
the weak-coupling regime. Under the rotation of the co- 
ordinate system, Ak changes its sign, if the Cooper pair- 
ing appears in the d-channel. The Cooper pairing might 
also take place with finite center-of-mass momentum, if 
electrons are spin polarized . 

On the other hand, the symmetry of the order parame- 
ter could be different from the 'internal' symmetry of the 
pair wave-function, and from the symmetry of a single- 
particle excitation gap in the strong-coupling regime |3| ■ 
Real-space pairs might have an unconventional symme- 
try due to a specific symmetry of the pairing potential as 
in the case of the Cooper pairs The d-wave symme- 
try of the ground state could be also due to a topological 
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degeneracy of inter-site pairs on a square lattice jsO], as 
proposed in Ref. [sJl- But in any case the ground state 
and DOS are homogeneous, if pairs are condensed with 
K = 0. However, if a pair band dispersion has its minima 
at finite K in the center-of-mass Brillouin zone, the Bose 
condensate is inhomogeneous. In particular, the center- 
of-mass bipolaron energy bands could have their min- 
ima at the Brillouin zone boundaries at K = (tt, 0) and 
three other equivalent momenta |32j (here and further I 
take the lattice constant a = 1, and fi = 1). These four 
states are degenerate, so that the condensate wave func- 
tion ■(/'(m) in the real (Wannier) space, m — {mx,my), is 
their superposition, 

^(m)= &Ke-'^^"", (1) 

K=(±7r,0),(0,±Tr) 

where 6k = iy/n2/'2 are c-numbers, and nc{T) is the 
atomic density of the Bose-condensate. The superposi- 
tion, Eq.(l), respects the time-reversal and parity sym- 
metries, if 

il^im) — ^/jTc [cos(7rm2:) ± cos(7rmy)] . (2) 

Two order parameters, Eq.(2), are physically identical 
because they are related by the translation transforma- 
tion. They have c?-wave symmetry changing sign in the 
real space, when the lattice is rotated by 7r/2. This sym- 
metry is entirely due to the pair-band energy dispersion 
with four minima at K ^ 0, rather than due a spe- 
cific pairing potential. It reveals itself as a checkerboard 
modulation of the hole density with two-dimensional pat- 
terns, oriented along the diagonals. From this insight one 
can expect a fundamental connection between stripes de- 
tected by different techniques [sM Is^ a nd the symmetry 
of the order parameter in cuprates l32| . 

Now let us take into account that in the supercon- 
ducting state (T < Tc) single-particle excitations interact 
with the pair condensate via the same short-range attrac- 
tive potential, which forms the pairs |23|. The Hamilto- 
nian describing the interaction of excitations with the 
pair Bose-condensate in the Wannier representation is 

H = - Y [t(m - n) + M'^m,n]4mCs„ 

s,m,n 

+ $^[A(m)c|^cx™ + H.c.], (3) 

m 

where s =|, | is the spin, t(m) and are hopping inte- 
grals and the chemical potential, respectively, and 
Csm create (annihilate) an electron or hole at site m, and 
A(m) c>c '(/'(m)- Applying equations of motion for the 
Heisenberg operators ct^ (t) and Cgmit), and the Bogoli- 
ubov transformation |35j 

5Tm(t) = 5][?/.(m)a,e-'^'^* + f:(m)/3te--*], (4) 



5im(t) = ^K(m)/3.e-"''* - <(m)ate-"*], (5) 



one obtains BdG equations describing the single-particle 
excitation spectrum, 

eu(m) = - ^[i(m - n)+/i(5m,n]w(n)-)- A(m)i;(m), (6) 



-ev{m.) = - ^[i(m - n) -I- ^(5m,n]w(n) -I- A(m)u(m), 

n 

. (7) 

where excitation quantum numbers v are omitted for 
transparency. These equations are supplemented by the 
sum rule X)i.[w;y(m)u*(n) -I- v^{m.)v*{n)] = (5m,n, which 
provides the Fermi statistics of single particle excitations 
a and f3. Different from the conventional BdG equations 
in the weak-coupling limit, there is virtually no feedback 
of single particle excitations on the off-diagonal poten- 
tial, A(m), in the strong-coupling regime. The number 
of these excitations is low at temperatures below Ap/fc^, 
so that the coherent potential A(m) is an external (rather 
than a self-consistent) field, solely determined by the pair 
Bose condensate p7j . 

While the analytical solution is not possible for any 
arbitrary off-diagonal interaction A(m), one can readily 
solve the infinite system of discrete equations (6,7) for 
a periodic A(m) with a period commensurate with the 
lattice constant, for example 

A(m) = AJe'''"="-e^''"«], (8) 

which corresponds to the pair condensate at K = (±7r, 0) 
and (0,±7r), Eq.(2), with a temperature dependent (co- 
herent) Ac oc ^JnJ^T). In this case the quasi-momentum 
k is the proper quantum number, = k, and the excita- 
tion wave-function is a superposition of plane waves, 

u,(m) = Wke^*'™ + uke^**'~^^"", (9) 

Here g^; = (tt, 0), = (0, tt), and g = (tt, tt) are recipro- 
cal doubled lattice vectors. Substituting Eqs.(9) and (10) 
into Eqs.(6,7) one obtains four coupled algebraic equa- 
tions, 

CkWk = CkWk - Ac(i'k - i^k), (11) 

EkMk = '?k-gWk + Ac(wk - Wk), (12) 
-EkWk = 'Ck-g^.fk + Ac(uk - Uk), (13) 
-EkWk = Ck-g^Wk - Ac(wk - Uk), (14) 

where ^k = ^ X^n ^(n)e*'' " — /i. The determinant of the 
system (11-14) yields the following equation for the en- 
ergy spectrum e: 

(e - ^k)(e - Ck-g)(e + Ck-gJ(e + ^k-gj 
= A2(2e-t-a-g. +ek-gj(2e-ek-ek-g). (15) 

Two positive roots for e describe the single-particle exci- 
tation spectrum. Their calculation is rather cumbersome. 
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FIG. 1: Single-particle excitation energy spectrum ( arb. 
units) along the diagonal direction of the two-dimensional 
Brillouin zone 



but not in the extreme strong-coupling limit, where the 
pair binding energy 2Ap is large compared with A^, and 
with the single-particle bandwidth w. The chemical po- 
tential in this limit is pinned below a single-particle band 
edge, so that fi = — (Ap -I- w/2) is negative, and its mag- 
nitude is large compared with Ac. Then the right hand 
side in Eq.(15) is a perturbation, and the spectrum is 



eik 



e2k 
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(16) 
(17) 



Its dispersion along the diagonal direction is shown in 
Fig.l in the nearest neighbor approximation for the hop- 
ping integrals on a square lattice. 

If a metallic tip is placed at the point m above the 
surface of a sample, the STM current I{V, m) creates an 
electron (or hole) at this point. Applying the Fermi-Dirac 
golden rule and the Bogoliubov transformation, Eq.(4,5), 
and assuming that the temperature is much lower than 
Ap/Zcs one readily obtains the tunnelling conductance 



cr(F, m) 



dI{V, m) 
W 



oc 



\u,{m)\''5[eV -e. 



(18) 



which is a local excitation DOS. The solution Eq.(9) leads 
to a spatially modulated conductance, 

(t{V, m) = (Jreg{V) + l7raod{V) COs(7rm^ -f TTinrLy) . (19) 

The smooth (regular) contribution is 



Co > [U^ 

k,r=l,2 



+ ui^)5{eV -er^,), (20) 



and the amplitude of the modulated contribution is 

iyr,iod{V) = 2(Ta ^ Wrk'"rk'5(ey - e^k), (21) 



where ctq is a constant. Conductance modulations reveal 
a checkerboard pattern, as the Bose condensate itself, 
Eq.(2), 



reg 



(22) 



where 



A = 2 ^ [uikUik<5(ey - eik) + U2kW2k'5(ey - e2k)] / 

k 

E [Kk + ^lidK^y - eik) + ("2k + "2k)'5(e^ - e2k)] 



is the amplitude of modulations depending on the voltage 
V and temperature. An analytical result is obtained in 
the strong-coupling limit with the excitation spectrum 
given by Eqs. (16,17) for the voltage near the thresh- 
old, eV « Ap. In this case only states near bottoms of 
each excitation band, Fig.l, contribute to the integrals 
in Eq.(22), so that 



and 



"Ik 



M2k 



Ck — eik 
eik — ^k-g 



-wik 
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k-g 



e2k — Ck 



e2k - 

W2k 



-Wlk < Ulk, 



-U2k < U2k- 



(23) 



(24) 



Substituting these expressions into A, Eq.(22), yields in 
the lowest order of A^, 



A: 



2A| 



(25) 



The result Eg. (2 2) g enerally agrees with the STM ex- 
periments j24l. I25I l26l l36l ISTjl . where the spatial checker- 
board modulations of a were observed in a few cuprates. 
The period of the modulations was found cither com- 
mensurate or non-commensurate depending on a sample 
composition. In ovir model the period is determined by 
the center-of mass wave vectors K of the Bose-condensed 
preformed pairs. While the general case has to be solved 
numerically [33|, the perturbation result, Eq.(22) is qual- 
itatively applied for any K at least close to Tc, where the 



coherent gap is small, if one replaces cos(7rm^ 



k,r=l,2 



by cos{Kxmx -\- Kyiriy). The period of DOS modula- 
tions does not depend on the voltage in the perturbation 
regime, as observed [s^, but it could be voltage depen- 
dent well below Tc, where higher powers of Ac are im- 
portant. Different from any other scenario, proposed so 
far ^3 J the hole density, which is about twice of the con- 
densate density at low temperatures, is spatially mod- 
ulated with the period determined by the inverse wave 
vectors corresponding to the center-of-mass pair band- 
minima. This 'kinetic' interpretation of charge modula- 
tions in cuprates, originally proposed [s^l before STM 
results became available, is consistent with the inelas- 
tic neutron scattering, where incommensurate inelastic 
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peaks were observed only in the superconducting state 
|4l| . The vanishing at of the incommensurate peaks 
is inconsistent with any other stripe picture, where a 
characteristic distance needs to be observed in the nor- 
mal state as weh. In our model the checkerboard charge 
modulations should disappear above Tc, where the Bose- 
condensate evaporates and the coherent gap Ac{T) van- 
ishes, so that A = Om Eq.(22). While some STM studies 
|39| report incommensurate DOS modulations somewhat 
above T^, they might be due to extrinsic inhomogeneities. 
In particular, preformed pairs in the surface layer could 
bose-condense at higher temperatures compared with the 
bulk Tc- Our model is microscopically derived using the 
strong-couphng (bipolaron) extension of the BCS theory 
Pi . If the electron-phonon interaction is strong, such that 
the BCS coupling constant A > 1, electrons form bipo- 
larons above Tc, which are Bose condensed below Tc. The 
polaron bandwidth is exponentially reduced, which ex- 
plains a low estimate of the Fermi energy using the exper- 
imental London penetration depth in cuprates jl^ . Evi- 
dence for an exceptionally strong electron-phonon inter- 



action in high-temperature superconductors is now over- 
whelming (see, for example, E3)- Y^*' generally, 
the model describes charge modulations due to the Bose 
condensation with non-zero center-of-mass momenta of 
preformed pairs formed by any pairing interaction. 

In conclusion, I solved BdG equations with the periodic 
off-diagonal potential caused by the Bose condensation of 
preformed pairs with non-zero center-of-mass momenta, 
and found the checkerboard modulations of the single- 
particle DOS similar to those observed in cuprates. The 
main assumption that the ground state of superconduct- 
ing cuprates is the Bose-Einstein condensate of preformed 
pairs, is supported by a growing number of other exper- 
iments 8]. The model links charge orders, pairing, and 
pseudo-gaps as manifestations of a strong attractive in- 
teraction in narrow bands. 
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